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Abstract 

Systems with very long-range interactions (that decay at large distances like 
U(r) ~ r~ l with I < d where d is the space dimensionality) are difficult to 
study by conventional statistical mechanics. Examples of these systems are 
gravitational and charged (non-electroneutral). In this work we propose two 
alternative methodologies to avoid these difficulties and capture some of the 
properties of the original potential. The first one consists in expressing the 
original potential in terms of a finite sum of hard-core Yukawa potentials. 
In the second one, the potential is rewritten as a damped potential, using a 
damping function with a parameter that controls the range of the interaction. 
These new potentials with finite ranges, which mimic the original one, can 
now be treated by conventional statistical mechanics methods. 
Keywords: Long-range interactions, discrete perturbation theory, 
first-order mean spherical approximation 
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1. Introduction 

Description of systems interacting via the so-called long-range interac- 
tions (LRI) is an important statistical mechanics problem. These systems are 
found from very small to very large scales, for instance, in astrophysics 
3|, plasma physics J4, 5|, hidrodynamics Jg|, atomic physics J?J and, nuclear 
physics Jf|. 

In order to use a precise definition of LRI between a pair of particles, 
that are a distance r apart, we consider the following: when the interaction 
potential between particles decays at long distances like l/r l in a space of d 
dimensions, the interaction can be considered to be long-range if I < d. 
This definition is a consequence of considering the energy e of a given particle 
located at the center of a sphere of radius R with a homogeneous particle 
distribution in ^-dimensions. In order to exclude the divergence that appears 
at very short distances, the energy of the neighboring particles located inside 
a sphere of radius 5 is neglected, e is given by, 

e = [ R Pl d d r = pB Q d [ R r^- l dr 
Js r Js 

d — a 

where p is the generic particle density, B is a coupling constant which guar- 
antees the correct energy dimensions, and is the angular volume in the 
d-dimensional space. When R is increased, e remains finite only when / > d; 
such cases are the usual short-range interactions. The opposite corresponds 
to I < d, where energy diverges for an increasing volume; these are long-range 
interactions. Examples of different long-range potentials are shown in Figure 
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1. Notice that this particular definition could be different in the context of 
fluids theory. 

In statistical mechanics, most of the effort to obtain the equilibrium and 
non-equilibrium thermodynamic properties, has been concentrated on sys- 
tems with short-range interactions. One of the main features of LRI systems 
is that their total energy, under the pairwise additive approximation, is non- 
extensive, and as a consequence, is also non-additive |9l-|l4|. Therefore, the 
connection between Boltzmann-Gibbs (BG) statistical mechanics and clas- 
sical thermodynamics is not straightforward, since the latter assumes that 
energy is an additive quantity 15j |. To our knowledge, there is not a thermo- 
dynamic formalism (independent of a statistical mechanics approach) that 
allows this connection; however, it is possible to start from a non-extensive 
statistical mechanics and to obtain a non-extensive thermodynamic formal- 
ism. 

Perhaps a non-extensive version of statistical mechanics could be a more 
natural theoretical frame to study LRI. A few proposals for BG statistical 



mechanics generalizations have been given 



tew p 



17] , however none of them 



are unanimously accepted. Besides, the application of these generalizations 
to long-range potentials has been scarce. Another approach is to make ad- 
justments to the BG formalism to study these systems J^, [nj, 18 |. 

In order to avoid the difficulties to treat LRI mentioned above, in this 
work we present a first naive approach, but general, in the sense that it can 
be applied to a great variety of long-range potentials in the frame of BG 
statistical mechanics. This methodology consists in rewritting a long-range 
potential as a short- range one, being the latter similar to the long-range 
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potential in its graphical representation. We expect that this short-range 
potential recovers some features of the original one, however we know that 
this path leads to a classical thermodynamics frame and we do not know if 
real systems with LRI are well represented by this thermodynamics. 
More specifically our approach consists in expressing a given long-range po- 
tential as a) a finite sum of Yukawa potentials; b) a product of the potential 
with a damping function which depends on a parameter, that under a certain 
limit, the original potential is recovered. 

We have selected the Discrete Perturbation Theory (DPT) and the First- 



order Mean Spherical Approximation 



20| (FMSA) to study these potentials. 



These theories have been successfully applied in the context of fluids and 
more recently in the soft matter field, and can be applied to a great variety 
of potentials. 

As an illustrative example, we have choosed the gravitational potential 
due to the interaction between two identical spherical rotating bodies (ETS 
potential), obtained by Escamilla et al. 2jJ. More interesting models in the 
context of molecular liquids, could be, for instance, the Coulomb interactions. 

This work is organized as follows. In section [2j we give a brief descrip- 
tion of the ETS potential and of the hard-core multi- Yukawa (HCMY) and 
damped potential approaches. In section [3] we present internal energies, pres- 
sures, and vapor-liquid phase diagrams for the approximated potential. Fi- 
nally, in section H] we give the main conclusions of this work. 
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2. Theory 

2.1. The ETS potential 

In the context of general relativity, within the weak-field limit methodol- 
ogy, an angular averaged potential due to the interaction between two iden- 
tical spherical rotating bodies was proposed 2JJ. This interaction potential 



for hard-core spheres is given by U*(x) = U(x)/ | e m j n |, where: 



oo if x < 1 

U*(x) = | (2) 

1 1 ? »■> arctan (a* Ix) if x > 1; 

with a* = J/Mca, M is the mass, J is the angular momentum, c is the 
speed of light in vacuum, a is the diameter of particles, x = r/a and e m j n 
is the potential evaluated at x = 1. This potential is purely attractive, 
nondivergent at short distances (for x — > , 1/a* arctan(a*/x) — > vr/2), and 
keeps its long-range nature satisfying the condition given by (OQ). Specific 
angular momentum a* is a parameter which modulates the intensity of the 
interaction. The long-range behavior is the same for any finite value of a*. 
For instance, in the limit of a* — > the ETS potential goes to conventional 
— 1/x Newtonian gravitational interaction. In Figure 2, the ETS and —1/x 
potentials are shown; it can be noticed that the long-range behavior is the 
same for both of them. 

To avoid the difficulties to evaluate thermodynamic properties for this 
long-range potential, we propose to rewrite the non hard-core potential part 
as: 

1) a finite sum of m Yukawa potentials, 
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= e* ; multi-Yukawa approach, (3) 

i=l 

with the energy and range parameters e* = 6j/ | e m j n | and «*, respectively, 
and 

2) a damped potential, which consists in the product of the original potential 
and a damping function f{^j,x), 

U(x) 

$b(7, x) = /(7, x)-. r\ damped approach, (4) 

I tmin I 

where 7 is the damping parameter that can be selected in order to guarantee 
that the approximated potential mimics the original one. 

2.2. First-order mean spherical aproximation 

?he first-order mean spherical approximation was developed by Tang et 



al., 22] as an improvement of the mean spherical aproximation (MSA) [23j|. 

The solution of the Ornstein-Zernike integral equation under MSA makes 
it possible to find analytical thermodynamic and structure expressions, which 
otherwise would require time-consuming numerical work. Despite these ad- 



vantages, MSA may lose its solution in unstable regions |24j. An improve- 
ment of this theory is the first-order mean spherical approximation. 
FMSA solves analytically the radial distribution function (RDF) to first or- 



der in terms of inverse temperature. So 



and always exist in unstable regions 



24 1 



utions obtained are explicit, simpler 



A successful application of FMSA theory was done by Tang et al. 



22] 



to the Yukawa and HCMY potentials. A finite sum of Yukawa potentials 



(3 



can mimic other well-known potentials, like Lennard- Jones potential 20] or 



sticky hard spheres 



25] ; however, its efficacy for LRI like the ETS potential 



has never been tested to our knowledge. In this work we will approximate 
the ETS potential using HCMY potential to express it. 

The hard-core Yukawa potential with multiple tails m, for hard-core par- 
ticles is given by, 



MY 



OO 



if x < a 



(5) 



Z^i=l H x 11 x — °i 

For a system of N particles confined in a volume V at temperature T 
interacting through a HCMY potential ([5]), the reduced Helmholtz free energy 
(a = A/NksT^kB is the Boltzmann's constant), within the FMSA can be 



expressed as 



20] 



a = a idea i(r), T) + a hs {rf) + 01(77, T, k*, e*) + a 2 ((r/, T, «*, e*); (6) 

where a i( i ea i(ri,T) and 0^(77) are the ideal gas and hard-sphere contributions 
to free energy, rj = np*/6 is the packing fraction and p* = pa 3 is the reduced 
density. The Carnahan and Starling EOS 1 for 0/^(77) is used, 



4?7 — 3rf 



(7) 



The first and second order perturbation terms, ai(r], T, k* , e*) and 02(77, T, ft* , e*) 
are defined as 



e*L(n* 



rx,o- 12^ E (1 _,)2g o( ;* )K * 



:2' 



(8) 
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e*e* 

* 3 



with /3 = 1/ksT, n = 1, . . . ,m. The other functions that appear in 
are defined by 

Qo{t) ~ — — ' 

S(t) = (1 - r/) 2 t 3 + 677(1 - r])t 2 + ISrft - 12^(1 + 277), 



(9) 
and 



L(t) 



1 + 



(10) 



Other thermodynamic properties can be derived from ([6]). 



Discrete perturbation theory 
The discrete perturbation theory was developed by Benavides and Gil- 



Villegas 



19| and has been successfully applied to different discrete and con 



tinuous potentials 



27H29j . In particular Torres- Arenas et al. 30(, have found 



that this approach accurately describes the thermodynamics of a hard-core 
attractive Yukawa potential with k* = 0.1, which represents a long-range 
potential in the context of molecular liquids; however, it is a short-range in- 
teraction in the sense of ([[]) . In this work, DPT will be used for even smaller 
values of k*. 

Analogous to FMSA, DPT provides a recipe for the Helmholtz free-energy 
of an arbitrary radial potential, re-expressing it as a sum of square-well 
and/or square-shoulder potentials. For a system of N particles, contained in 
a volume V at temperature T, the reduced free Helmholtz energy is given 

by, 



a = a h , 



i=l 



af (77, \ iy ei) - af(r), A;_i, q) 



1=1 



af(?7, \ h 6i) - af (77, Ai_i, e 4 ) 



(11) 



where n is the total number of steps, 6j is an energy parameter that can be 
positive or negative and defines the height/depth of the step, and Aj — Aj_i 
is the width of i-th step, af (rj, X n , e$) and af (77, A n , 64) are the first-order and 
second-order perturbation contributions for a square-well (e, > 0) or square- 
shoulder (e, < 0). The evaluation of ( ITT]) can be simplified since square- well 
(SW) and square-shoulder (SS) potentials differ only in the sign of q and the 
following relations are satisfied: 

af S (v, e i) = - a i W (v, K 

af(r / ,A 4 ,e J ) = af iy (r / ,A l ,e J ). (12) 

The problem is now reduced to calculate the free energy perturbation terms 
of the square-well potential. 

In order to apply the DPT to a continous potential, a discrete version 



of the potential is required |30| . We used an approximation of the desired 
continous potential (either the HCMY or the damped potential approach) in 
terms of step functions, choosing the middle point on each step to evaluate 
the potential. 

For monotonic potentials without a defined range, DPT requires a cutoff; 
this cutoff can be calculated by the condition |?7(A C )| = 10 -6 , where A c 
is the cutoff distance. This condition was successfully applied to Yukawa 
potential 30 1. 



In the DPT approach, the perturbation terms af w (r], Aj, e,) and af w (rj, Aj, e, 
have been calculated from two theoretical equations: for 1.1 < A < 3.0, we 
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used Espmdola 



expressions 



et al., 3l| and, for longer SW ranges, Benavides and del Rio 



32] . DPT requires the SW equation of state for all values of A; 
however, the EOS proposed by Espmdola et al., have problems for A < 1.1. 
Therefore, in this work the number of discretizations, n, is constrained to pre- 



vent an evaluation of 



al w in 1.0 < A < 1.1 



me perturbation contributions for free energy af w and 



30[, n is determined by n = (A c — Aq)/0.1 = 10(A C — Aq). 



3. Results 

3.1. Multi-Yukawa approach 

The first proposal to study the thermodynamic properties of a long-range 
potential is to rewrite it as a finite sum of Yukawa potentials. This goal can 
be achieved usin g th e Levenberg-Marquardt method to perform a non-linear 



least squares fit 



33|. 



For the ETS potential, we have found that three Yukawa potentials are 
enough to describe it (see Figure 3). Notice that for both potentials, the 
long-range behavior is essentially the same and differences can be seen only 
at short distances. We expect that these short-range dicrepancies will be 
irrelevant for this type of potentials. The fitting was performed for the par- 
ticular case of a* = 10. 

The ETS potential in the HCMY approach is given by ([5]) for m = 3. 
Corresponding parameters e* and k* are presented in Table [TJ 

As can be seen, k* values are small, corresponding to long-range molec- 
ular potentials; nevertheless, they are short-ranged if we consider definition 
given in ([I]). Therefore, HCMY potential can be treated by the conventional 
methods of BG Statistical Mechanics. 
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Table 1: Strength and range parameters for HCMY approach to the ETS potential. 

i e* k* 



1 5.6590 0.1500 

2 -2341.2060 3.6744xl0~ 4 

3 2334.5470 3.6880 xlO" 4 



In the frame of a linear theory, as is the case of FMSA, the problem is reduced 
to the knowledge of an EOS for the hard-core attractive/repulsive Yukawa 
potential of variable k* . Yukawa EOS have not been tested for such small 
k* values. Results obtained in this work can also be used as a test of the 
reliability of FMSA and DPT in a domain of «j « 10 -4 , far from the n* 
values commonly studied by liquid state community 
The reduced units that will be used are given by, 



22 



30 



34|. 



k B T 



Pa' 



| ^min | | ^min \ KB-*- 

where T*, P*, fi*, and p* are the reduced temperature, pressure, chemical 
potential, and density, respectively. 

To estimate the accuracy of FMSA, compared with DPT for a hard-core 
attractive Yukawa potential (HCAY) with small values of the range param- 
eter, we considered the case k* = 0.1. For this value, DPT gives very good 



predictions 



30|. To our knowledge, it is the smallest value considered in sim- 



ulation studies 



35|. 



The results for a vapor-liquid phase diagram with FMSA, DPT and simu- 
lation are presented in Figure 4. Both theories are in good agreement with 
simulation data. 
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With the performance test of FMSA and DPT done, we constructed 
vapor-liquid phase diagrams of the HCMY approach to the ETS potential 
for a* = 10. In Figure 5, vapor-liquid phase diagrams for both theories are 
presented. 

In the absence of simulation data for the thermodynamic properties of the 
ETS potential, we present only a theoretical prediction of pressure and excess 
internal energy as a function of density, for several isotherms. 

The cutoff (A c ) for the potentials used within DPT has been selected as 
|^(Ac)| = 10~ 8 that guarantees that the thermodynamic properties do not 
change significatively if the cutoff is increased. 

In Figure 6 and in Figure 7 we present reduced pressure and excess in- 
ternal energy as a functions of density, for the approximated predictions of 
FMSA and DPT theories of HCMY. In Figure 6 the selected isotherms were 
chosen around the critical temperature. For the case a* = 10, coexistence 
isotherms (showing a van der Waals loop) appear at very hig h temperatures; 



30 



35|. 



which is a typical behavior for long-range potentials 
The corresponding excess internal energies shown in Figure 7, present the 
same ordinary functionality as the one exhibited by a short-range potential. 

3.2. Damped potential approach 

In order to illustrate the methodology given in (jlj) for the ETS potential, 
we have selected the damping function, /(7, x) = exp[— 7(2; — 1)]. Therefore 
the damped potential for x > 1, can be written as: 



arctan(a *) \ x 



12 



where we considered the case a* = 10. The parameter 7 = 1.9603 x 10 
was estimated by solving, 



This equation guarantees that, in this approach, long-range first-order per- 
turbation term (g(x) ~ 1, g(x) being the radial distribution function) is 
equivalent to the one in the HCMY approach. In this way we avoided the di- 
vergence of this term, if the ETS potential is used. It is important to remark 
that the previous condition is not the only possible one; for instance, another 
potential that mimics the ETS potential expressed in a different basis, could 
have been chosen instead of &my(%)- 

With this 7 selection, in Figure 8 we show that damped and the ETS po- 
tentials are very similar. This new potential will be studied using DPT. Since 
there are no available simulation data for the ETS potential, in the following 
figures we make a comparison between the two methodologies presented in 
this work. 

In Figure 9 the vapor-liquid phase diagrams of damped potential and 
HCMY approximations to the ETS potential are shown. As can be seen 
they give very similar predictions. 

In Figure 10 and Figure 11, reduced pressures and excess internal energies 
as a function of density for damped potential and HCMY approaches are 
presented. It can be noticed that both methodologies are equivalent. 




(15) 
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4. Conclusions 

In this work two different methodologies to study long-range potentials 
fluids from the point of view of BG statistical mechanics have been pre- 
sented. These methods, HCMY and damped potential, are equivalent when 
used to predict the thermodynamic properties of a selected long-range poten- 
tial; however, simulation data is required to check their accuracy. A similar 
treatment can be applied to other long-range potentials. 
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Figure captions 

Figure 1. Examples of different (repulsive and attractive) long-range 
potentials, that at long distances decay as 1/r. 

Figure 2. The ETS potential with a* = 1 (solid line) compared with 
— 1/x potential (dashed line). Long-range behavior is the same for both po- 
tentials. 



Figure 3. Comparison between the ETS potential for a* = 10 (solid 
line) and HCMY potential (dashed line) obtained with the non-linear least 
squares fitting, with k* and e* given in Table [TJ 



Figure 4. Vapor-liquid phase diagram for a HCAY for k* — 0.1. At 
this value of k*, FMSA and DPT predictions overlap. Diamonds are MC 



simulation data from Caillol et al. 



35|. 



Figure 5. Vapor-liquid phase diagram for HCMY approach to the ETS 
potential. Solid line represents DPT prediction, and squares represent FMSA 
prediction for HCMY approach with e* and «*, given in Table [TJ 



Figure 6. The ETS potential pressure-density plot within HCMY ap- 
proach for T* = 1.2 x 10 8 (bottom), T* = 1.7 x 10 8 (middle), both subcritical 
temperatures, and a supercritical, T* = 3.0 x 10 8 (top) for DPT (solid line) 
and FMSA (dashed line). 
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Figure 7. The ETS potential reduced excess internal energy within 
HCMY approach for DPT (solid line) and FMSA (dashed line) at two su- 
percritical temperatures T* = 2 x 10 8 (bottom) and T* = 3 x 10 8 (top). 

Figure 8. Comparison between the ETS potential for a* = 10 (solid 
line) and the damped potential with damping parameter 7 = 1.9603 x 10~ 4 
(dashed line). 

Figure 9. Vapor-liquid phase diagram prediction within DPT of damped 
approach considering 7 = 1.9603 x 10 -4 and HCMY with e* and k* given in 
Table [TJ the prediction under DPT for both approaches overlap. 

Figure 10. Reduced pressure as a function of reduced density using DPT 
for damped potential (dashed line) with damping parameter 7 = 1.9603 x 
10 -4 , and HCMY potential (solid line). Two temperatures are depicted, su- 
percritical T* = 3.0 x 10 8 (top) and subcritical T* = 1.5 x 10 8 (bottom). 

Figure 11. Reduced excess internal energy as a function of reduced 
density using DPT for damped potential approach with a damping parameter 
7 = 1.9603 x 10~ 4 (dashed line) and HCMY (solid line). Two supercritical 
isotherms are presented, T* = 3.0 x 10 8 (bottom) and T* = 5.0 x 10 8 (top). 
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Figure 5. 
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Figure 9. 
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Figure 11. 
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